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1.  Introduction 

Several  interesting  mathematical  problems  concerned  with  partition  into 
classes,  and  surface  covering  are  suggested  by  the  physical  mechanisms  and 
geometry  of  antibody  attachment  to  virus  particles.  This  paper  outlines  some 
of  the  recent  work  done  in  this  area  by  Yassky  [7],  myself  [2],  [3],  Moran  and 
Fazekas  de  St.  Groth  [5]  and  Gilbert  [4],  adds  a  model  and  other  extensions 
of  my  own,  and  concludes  with  a  suggestion  for  further  investigations.  I  shall 
endeavor  throughout  to  hold  the  mathematical  argument  at  a  simple  level,  and 
emphasize  the  model  building  aspect  of  the  work,  in  the  hope  that  virologists 
may  be  tempted  to  use  and  perhaps  verify  experimentally  some  of  the  models 
put  forward. 

Let  us  consider  at  any  time  t  ^  0,  a  nutrient  medium  (either  in  the  laboratory, 
or  within  a  living  animal)  in  which  there  exist  a  fixed  number  N  of  particles  of 
a  virus  V;  suppose  that  at  time  t  =  0,  there  are  x0  >  N  antibodies  released  into 
this  medium.  We  may  expect  the  antibodies  to  attach  themselves  progressively 
in  some  random  fashion  to  the  viruses,  both  types  of  particles  being  subject  to 
Brownian  motion.  If  each  virus  particle  permits  a  maximum  of  s  attachments, 
then  at  any  time  t  ^  0,  the  N  virus  will  be  divided  into  s  +  1  classes  consisting 
of  n0(t),  ■  •  •  ,  ns(t)  particles  with,  respectively,  0,  1,  •  •  •  ,  s  antibodies 

attached  to  them;  there  will  remain  x(t)  =  x0  —  X)Ui  ini(t)  unattached  antibody 
particles.  The  constitute  a  class  partition  of  the  virus  particles,  which 

varies  in  time  t.  We  may,  for  simplicity,  in  some  cases,  approximate  the  integer 
valued  random  variables  0  ^  w,(£ )  ^  N,  and  0  ^  x(t)  ^  x0  by  analogous  func¬ 
tions  differentiable  in  t;  then,  as  we  shall  see,  a  deterministic  approximation  to 
the  random  evolution  of  the  {n»(0}  and  x(t )  can  be  found.  It  is  also  possible  to 
obtain  a  stochastic  approximation  to  the  integer  valued  using  the  previous 

deterministic  approximation  for  x(t). 

While  such  results  may  indicate  the  number  of  antibody  attachments  to  the 
virus,  they  do  not  alone  provide  adequate  information  as  to  its  loss  of  infectivity. 
Two  cases  arise,  however,  in  which  fuller  informal  ion  can  be  obtained.  These  are 
the  cases  when 
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(1)  the  particle  of  virus  V  is  approximately  spherical,  as  it  is  for  influenza,  and 

(2)  the  virus  V  is  a  bacteriophage. 

In  case  1,  Moran  and  Fazekas  de  St.  Groth  [5]  have  shown  that  the  attach¬ 
ment  of  antibodies  to  a  single  approximately  spherical  virus  particle  may  be 
formulated  as  a  problem  of  geometrical  probability.  Consider  i  cylindrical  an¬ 
tibodies  in  the  shape  of  long  cigars  adhering  by  one  of  their  ends  to  a  particle 
of  spherical  virus;  each  antibody  when  standing  normally  to  the  virus  surface 
shields  a  circular  spherical  cap  (subtending  a  half  angle  a  at  the  center  of  the 
sphere)  on  it  from  contact  with  a  healthy  cell.  For  the  influenza  virus,  the  radius 
of  the  sphere  is  40  m/u,  while  the  antibodies  are  of  length  27  m/z  so  that  the 
shielded  area  subtends  a  half  angle  a  =  53.43°.  A  sufficiently  large  number  i  of 
such  antibodies  would  result  in  a  complete  covering  of  the  sphere  and  cause  loss 
of  infectivity  of  the  virus.  Moran  and  Fazekas  de  St.  Groth  have  obtained  the 
asymptotic  value  for  large  i  of  the  probability  P(i)  that  the  sphere  is  covered 
by  i  antibodies,  and  more  recently  Gilbert  [4]  has  found  general  bounds  for 
this  probability.  If  at  any  time  t  ^  0  we  also  know  the  partition  of  the  N  virus 
particles  into  the  classes  {n,}  of  particles  carrying  i  =  0,  1,  •  •  •  ,  s  antibodies, 
we  can  evaluate  the  probability  of  loss  of  infectivity  of  the  virus  at  time  t. 

In  case  2,  it  is  known  that  a  single  antibody  attachment  to  the  bacteriophage 
tail  causes  loss  of  infectivity.  This  means  that  of  the  s  possible  attachments,  a 
single  particular  one  will  suffice  to  prevent  infectivity.  Suppose  we  now  parti¬ 
tion  the  rii  viruses  with  1  ^  i  ^  s  —  1  antibody  attachments  into  two  classes 

(1.1)  0  ^  ^  N,  0  n0i  ^  N,  +  nQi  =  rii, 

where  the  first  suffix  indicates  attachments  to  the  tail  and  the  second  to  any 
other  position  on  the  virus.  Clearly,  ni,,_i  =  n8,  n00  =  n0.  We  shall  show  that  the 
{fti.t-i},  {ftoi},  with  i  =  1,  •  •  •  ,  s  —  1  can  be  approximated  by  deterministic 
values,  and  also  obtained  stochastically  using  the  previous  deterministic  ap¬ 
proximation  to  x(t).  Loss  of  infectivity  in  the  deterministic  case  results  when  all 
virus  particles  have  acquired  tail  antibody  attachments;  in  the  stochastic  case 
the  probability 

(1.2)  P{n co  =  0,  ttoi  =  0,  •  •  •  ,  n0,s-i  =  0} 

will  give  a  measure  of  the  noninfectivity.  We  now  proceed  to  consider  these 
cases  in  detail. 

2.  The  attachment  pattern  of  antibody  to  virus 

Suppose  we  consider  first  a  deterministic  approximation  to  the  attachment  of 
antibodies  to  virus  particles  with  s  possible  emplacements.  For  time  t  ^  0,  let 
the  differentiable  functions  n0(t),  •  •  •  ,  n8(t)  approximate  the  number  of  particles 
with  0,  1,  •  •  •  ,  s  virus  attachments,  where  wt(<)  =  N ;  also  let  x(t)  =  x0  — 
be  the  antibodies  remaining  unattached  at  time  t.  Then  if  \t,  with 
i  =  1,  •  •  •  ,  s  (X,  =  0),  represents  the  attachment  rate  of  a  further  antibody 
to  a  virus  particle  already  carrying  i  of  these,  it  is  readily  seen  that 
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driQ 

dt 


—  \oTIqX, 


(2.1) 


drij 

dt 


=  (Xt_in4_i  —  X»n,)x, 


dx  *  - 1 

-tt  =  -xH 
at  t =o 


represent  the  state  equations  for  the  {n»(0}>  x(t).  The  initial  conditions  are 
(2.2)  Wo(0)  =  N,  n*( 0)  =0,  i  —  1,  •  •  •  ,  s,  rr(0)  =  x0. 

Writing  (2.1)  in  matrix  form,  we  obtain  following  the  transformation  p(t)  = 
x(t)  dr,  that 


(2.3) 


dn  _ 

* -  -Ln’ 


where  n'  =  (n0,  nh  •  •  •  ,  nt)  denotes  the  row  vector  of  the  {rii(t)}  and 

Xo 

—  Xo  Xi 


(2.4) 


L  = 


—  X„_2  Xg_x 

Xg_!  Xg  J 


If  L  is  written  in  the  canonical  form 


(2.5)  L  =  A_1AA, 

where  A  is  the  diagonal  matrix  of  eigenvalues  X0,  •  •  •  ,  X,,  assumed  distinct  for 
simplicity,  then  it  is  readily  shown  that  the  solution  to  (2.3)  is 

(2.6)  n(0  =  e_Lp(i)n(0)  =  A_1e_Ap(t)An(0) 

where  n'(0)  is  the  row  vector  (N,  0,  •  •  •  ,0).  The  equation  for  dx/dt  in  (2.1) 
becomes 


(2.7)  ^  - - xX'A-1e-Ap(t)An(0), 

where  X'  =  (Xo,  •  •  •  ,  X,),  and  this  may  be  rewritten  as 

(2-8)  S  -  ~di  VA’Ie ■i,<‘)An(0)  -  -I {.?„ 

the  Ci  being  constant  coefficients  of  e_x,p(e).  The  explicit  solution  of  this  second 
order  differential  equation  would  provide  the  complete  values  for  the  elements 
of  the  vector  n (t) ;  however,  such  a  solution  cannot  in  general  be  found. 

A  particular  case  in  which  (2.8)  can  be  solved  was  considered  in  a  slightly 
different  context  by  Yassky  [7].  Yassky  postulated  a  linear  relationship  for  the 
attachment  parameters  of  the  type 

(2.9)  X,-  =  a(s  —  i ), 


a  >  0, 
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and  found  that  the  system  then  simplified  considerably.  From  (2.1)  and  the 
relation  Xq  =  x  +  Y\~i  he  obtained  for  x(t)  the  differential  equation 

(2.10)  ~  +  x  Y  <*(s  —  i)ni  =  -jr  +  a{x  —  N(m  —  s)}x  =  0, 

ul  %  s  0  Wt 


where  m  =  x0/N  >  1  denotes  the  multiplicity  of  antibodies. 
The  solution  to  (2.10)  is 


(2.11) 


x(t)  = 


xp(s  —  m) 
—  m 


where  p  =  Na(s 
found  to  be 

no(0 


(2.12) 


««(0 


ns(t ) 


m),  and  from  it,  the  solutions  for  the  n*(f)  are  directly 
Ar{e_M<  +  s(s  —  m)_1(l  ~  e~M‘)}-*> 

(»)"*  {(^)  _1}’  *  =  i,  —  i, 

JV-’e 

1=0 


If  we  assume  the  attachment  of  antibodies  to  occur  as  a  Markov  process,  it 
is  simple  to  obtain  the  forward  Kolmogorov  equations  for  the  probabilities 
P(n0,  •  •  •  ,  ns]  x)  t)  that  at  time  t,  the  viruses  are  partitioned  into  classes  {«,} 
and  there  are  x  unattached  antibodies.  Although  such  equations  (of  birth  process 
type)  can  be  solved  in  principle,  they  prove  to  be  rather  intractable  in  practice, 
and  a  simplification  is  helpful.  This  consists  in  considering  the  reduced  stochastic 
process  for  which  x(t)  is  a  deterministic  differentiable  function  of  t,  while  the 
{nt}  are  stochastic  variables. 

Let  {Xo,  •  •  •  ,  X,}  be  the  set  of  attachment  parameters  such  that  Xi  >  0,  for 
i  =  0,  •  •  •  ,  s  —  1,  but  Xg  =  0.  The  probability  of  attachment  in  time  8t  of  an 
antibody  to  a  virus  already  carrying  i  phages  is  assumed  to  be 

(2.13)  \iUix8t  +  o(8t), 

where  rii  ^  0  is  the  random  number  of  bacteria  having  i  attached  phages, 
i  =  0,  ••*,«,  and  x{t)  =  dp/dt  is  the  deterministic  solution  of  equation  (2.8). 

Writing  Q  =  Q(n0,  nh  •  •  •  ,  n„;  t)  as  the  probability  that  at  time  t  ^  0  there 
are  no,  •  •  •  ,  ns  viruses  with  respectively  0,  •  •  •  ,  s  attached  antibodies,  we  find 
for  the  forward  Kolmogorov  equations 

(2.14)  ^  =  Y  X<n ixQ 

ot  1=0 

a-l 

+  Y  X*(w<  +  l)zQ(n0,  •  •  •  ,  n,-  +  1,  ni+ 1  —  1,  •  •  •  ,  n,;  t). 

i  =  0 


If  <p(u0,  •  •  *,  us;  t)  is  the  probability  generating  function  of  these  probabilities, 
then  equations  (2.14)  lead  to 


d<p 

~dt 


F  Xtx(wi+i  U{) 

t=0 


dUi 


(2.15) 
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This  is  a  particular  case  of  the  multivariate  Markov  process  originally  discussed 
by  Bartlett  [1]. 

I  was  able  to  show  [2]  that  for  a  general  nonincreasing  function  .r  =  x(t)  for 
unattached  phages,  the  partial  differential  equation  (2.15)  can  be  solved  to 
obtain  the  probability  generating  function  explicitly  as 

(2.16)  <p(uo,  •  •  •  ,ua;t)  =  ^Z)Wtal(i)j'  » 

with  probabilities 


(2.17) 


Q  = 


N\ 


n0!  * ' •  n„\ 


ad 


a 


n, 

s 


of  multinomial  form,  where  the  probabilities  a,-(<)  are  given  by 


(2.18) 

a,(0  =  Z  BojAjie  Xip(t), 

/- o 

i  =  0,  •  • 

•  ,  S, 

with 

Boo  —  1,  Boj  =  II  Xr(\r  Xy)  1, 

r  =  0 

j  =  1,  •  • 

•  ,  S, 

(2.19) 

A  a  =  1,  Aji  —  II  Xr_i(Xr  Xy)-1, 

r=j+ 1 

1 

•<S> 

o 

II 

l;*-  l,  •• 

■  ,  s, 

elements  of  the  matrices  B  =  A-1,  A,  respectively. 

The  expectations  for  this  process  are 

(2.20)  Eirii )  =  Ndi(t), 

and  the  variances  Var  (nd  and  covariances  Cov  (n,,  n3),  i  ^  j,  are  also  easily 
obtained.  If  the  X»  take  the  special  form  (2.9)  suggested  by  Yassky,  then 

(2.21)  p(t)  =  fo  x(r)dr  =  I  log  (S 

and  the  expectations  iVa,(<)  reduce  to  the  expressions  (2.12)  found  for  the  fully 
deterministic  case. 


3.  The  covering  of  spherical  virus  particles :  loss  of  infectivity 

The  problem  of  covering  a  spherical  virus  particle  by  a  sufficient  number  of 
cylindrical  antibodies  standing  normally  to  its  surface,  thus  preventing  virus 
contact  with  healthy  cells,  has  been  outlined  in  section  1.  We  saw  that  this  was 
reducible  to  the  geometrical  problem  of  covering  a  sphere  randomly  by  circular 
caps,  each  cap  subtending  a  half  angle  &  at  its  center. 

Moran  and  Fazekas  de  St.  Groth  have  pointed  out  in  their  paper  [5]  that  the 
problem  is  a  generalization  of  Stevens’  [6]  random  distribution  of  i  arcs  of 
length  xona  circle  of  unit  circumference,  for  which  the  asymptotic  probability 
of  coverage  for  large  i  is 
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(3.1)  P(i)  ~i(  1  -  xY~\ 

Using  extremely  ingenious  approximation  methods,  and  assuming  a  <  90° 
to  be  small  or  of  moderate  size  (as  in  the  case  where  a  =  53.43°  for  the  sphere 
of  radius  40  myu,  and  an  antibody  of  length  27  m m),  and  the  number  of  uncovered 
regions  to  follow  a  Poisson  distribution,  Moran  and  Fazekas  de  St.  Groth 
obtained  the  asymptotic  probability  of  coverage  for  large  i  as 

(3.2)  P( i)  ~  exp  -  Y2  [{2  [l  (1  +  eos<S)]'(l  +  £  tan*|  <s))_1  -  l]"‘- 

More  recently  Gilbert  [4]  has  derived  for  a  =  90°  the  exact  result 

(3.3)  P(i)  =  1  -  (i2  -  i  +  2)2-i,  i  ^  1, 

and  has  obtained  quite  generally,  for  any  angle  a  ^  90°  when  i  >  (sin2  |a)_1  + 
1,  the  bounds 

(3.4)  1  -  |  i(i  —  1)^1  —  sin2  ^  d'j  sin2  ^ 


^  P(i)  -  (l  -  sin2  ^ 


for  the  probability  P(i)  of  coverage  of  the  sphere. 

These  results  together  with  a  knowledge  of  the  {n,}  discussed  in  section  2, 
allow  us  to  consider  changes  in  the  loss  of  infectivity  with  time.  For  the  case 
where  the  (n,(J)}  are  deterministic,  we  might  for  convenience  count  their 
values  as 

n\  —  0  if  rii(t)  < 


(3-5) 


n\  =  j  if  j  -  ^  ^  rii(l)  <j  +  £  j  =  1,  •  •  ■  ,  N  -  1, 
«{  =  JVifJV  —  ^  fii(t)  S  N. 


In  this  case,  the  probability  of  loss  of  infectivity  Q/>(0  at  time  t  ^  0  will  be 
given  approximately  by 

(3.6)  QD{t)  =  fl  [P(i)]ni' 

t  =o 

where  n't  are  the  values  of  the  rq(£)  counted  as  above. 

Clearly,  coverage  of  the  virus  is  impossible  for  i  <  r,r  being  the  minimum 
number  of  antibodies  which  can  totally  cover  the  sphere.  Thus,  for  i  =  0,  •  •  •  , 
r  —  1,  P(i)  =  0,  when  —  n[  =  -  •  •  =  n'_i  =  0,  we  have  that 

(3.7)  QD(t )  =  n  [P(t)]< 

i  —  r 

These  results,  though  clearly  approximate,  will  give  some  indication  of  the 
progressive  loss  of  infectivity  of  the  virus  particles  in  time. 

For  the  stochastic  case,  the  probability  of  loss  of  infectivity  Qs(t )  is  given  by 
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(3.8)  Qs(t)  =  £  II  [F(*')]w,'Q(no,  •  •  •  ,ns;t) 

{«,•)  i  =  0 

=  E  [JP(0)]~  •  •  •  [*>(<>)]"•  „  ,  —  -  I  of  •  •  •  «?• 

=  [ooP(0)  +  •  •  •  +  aaP(s)]N 
=  [arP(r)  +  •  •  •  +  aaP(s)]N, 

where,  in  general,  the  a»(<)  are  those  given  in  (2.17)  and  the  Ndi(t)  reduce  to 
(2.12)  in  the  special  case  where  we  set  X,  =  (s  —  i)<x. 

We  illustrate  the  preceding  stochastic  process  by  means  of  an  elementary 
example.  Let  a  =  90°,  N  =  10,  x0  =  60 (m  =  6),  s  =  5,  a  =  1;  then  u  =  —10, 
and 

(3-9)  *(0  =  6  z^rnr.- 

It  follows  that  the  {«,-(<)}  are  given  by 

ao(<)  =  En0(t)/N  =  (6e10‘  -  5)~5, 

«i(0  =  Eni(t)/N  =  5(6e10*  -  6)(6e10'  -  5)~5, 

(3.10)  <n(0  =  En2(t)/N  =  10(6e10‘  -  6)2(6e10t  -  5)~5, 

o,(0  =  En3(t)/N  =  10(6e10t  -  6)3(6e10‘  -  5)~5, 

O4(0  =  Ern(t)/N  =  5(6e10<  -  6)4(6e10<  -  5)~5, 

0,(0  =  En,(t)/N  =  (6e10‘  -  6)5(6e10‘  -  5)~6. 

Since  r  =  4  in  this  case,  it  follows  from  (3.8)  that 


(3.11) 


<M0  =  5 


(6e10t  -  6) 4 

(6el°i  _  5)5 


P(  4)  + 


(6c1  °'  -  6) 5 

(6e10t  -  5) 5 


P(5) 


=  |  (3  -  &-■“) 


(6  -  6e-lo04 

(6  -  5<r10‘)5' 


This  provides  some  indication  of  the  dependence  on  time  of  the  loss  of  infec- 
tivity.  Two  remarks  seem  in  order.  First,  it  is  clear  that  since  for  &  90°  the 

results  given  by  Moran  and  Fazekas  de  St.  Groth  are  asymptotic  for  large 
values  of  i,  it  is  necessary  for  good  approximations  that  the  number  s  of  emplace¬ 
ments  be  large.  Secondly,  while  for  simplicity  in  the  model,  we  have  allowed 
random  attachment  of  the  antibody  in  any  position  on  the  spherical  virus 
surface,  there  are  in  fact  only  a  fixed  number  of  emplacements  with  specific 
positions  on  the  spherical  surface  at  which  the  antibody  may  adhere.  In  our 
example,  taking  s  =  5,  it  is  possible  as  t  — *  «  to  reach  the  limiting  probability 
Qs  =  5/16  of  noninfectivity.  This  would,  in  practice,  be  uselessly  small.  In  fact 
if  there  were  only  five  emplacements  on  the  virus  particle,  total  coverage  would 
occur  with  probability  1  with  5  attachments.  Thus,  while  the  proposed  model  is 
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not  entirely  unrealistic,  it  is  at  best  a  rough  approximation  to  the  true  structure 
of  the  process. 


4.  Antibody  attachment  to  the  tail  of  a  bacteriophage 

Let  us  now  suppose  that  the  virus  V  is  a  bacteriophage,  and  that  a  single 
antibody  attachment  to  its  tail  would  prevent  infectivity.  We  shall  for  simplic¬ 
ity  consider  the  case  where  the  general  attachment  parameter  X*  is  of  the  form 

(4.1)  X,  =  (s  —  i)a  i  =  0,  1,  •  •  •  ,  s 

as  in  (2.9),  though  the  subsequent  methods  apply  quite  generally  for  any  X,  ^  0. 

Considering  either  the  deterministic  or  the  stochastic  case,  we  note  that  if  tail 
attachments  are  not  distinguished  from  attachments  in  other  positions  of  the 
phage,  then  the  {wt-(/)}  of  (2.12)  or  the  probabilities  Q  of  (2.17)  will  fully  describe 
the  attachment  process. 

If,  however,  we  wish  to  distinguish  tail  attachments  from  others,  then  we 
must  concern  ourselves  with  the  classes 


{w00} , 

(4.2)  {woi}  ,  ,  t  1,  1 '  ■  ,  s  lj 

{Wl,s— l}  > 

of  bacteriophage  with  0,  i,  and  s  total  attachments,  respectively,  the  first  suffix 
position  indicating  a  tail  attachment.  We  see  that 

(4.3)  woo  =  wo,  w0i  +  Wi.v-i  =  w  i,  Wj  ,s_i  =  ns. 

Let  us  first  examine  the  deterministic  case.  Here,  the  attachment  parameters 
associated  with  w00,  noi,  i  =  1,  •  •  •  ,  s  —  1,  are  now  in  the  form: 


transition 

(0,  0)  ->  (0,  1)  or  (1,  0) 
(4.4)  (0,  i)  (0,  i  +  1) 


parameter 
X0  =  as, 


(s  -  1  -  i) 

(s  -  i) 


Xt 


=  a(s  —  1  —  i), 


(o,  i)  (i,  o 


X, 

(s  -  i) 


Thus,  we  may  write 


dwoo 

dt 


dn0 

dt 


—  asxnoo, 


dw0  . 

7/7  1 


a(s  —  f).TWo,-  1  —  iv(s  —  fi.TWo,-, 


dn0 

dt 


s 


—  1  =  Ctt'Wo,-2 


ax'Wos— i, 


(4.5) 


i  -  1 ,  •  •  •  ,  s  -  2, 
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or  in  matrix  form 

n0, 

where  p(t)  =  J*  x(t)  dr  as  before,  and  no  =  (n0o,  •  •  •  ,  n0*_i). 

Using  methods  similar  to  those  of  section  2,  the  solution  of  this  set  of  differen¬ 
tial  equations  is  easily  found  to  be 

(4.7)  n0i(t)  =  N  (^S  T  M  e-(«-*W(i  _  e~apy ,  i  =  0,  •  •  •  ,  s  -  1, 

where  p(t)  =  or1  log  {(s  —  me~tit)/(s  —  m )}  as  in  (2.21).  It  follows,  therefore, 
since 

(4.8)  nm(t)  =  n0(t)  =  N^er*1  -f  s(s  —  m)~1(  1  —  e_,i0}~s 
as  before,  that  the  numbers  of  phage  with  tail  attachments  are 

fli,i_i(0  =  »*(0  —  noi(t) 

(4’9)  =(i-l)"°{(^)"*  “1}' 


s-l 

8-1  /, 

a  r  a 

nn\  —  8_1 

wi,s_i  =  ns  =  N  —  Y.  n( 

II 

1 

M 

.)  n0  j  (■ 

4)  -q 

II 

o 

t=0  V 

*/  IV 

N  /  J 

Immunity  may  be  considered  to  exist  if  all  phage  have  a  tail  antibody  attach¬ 
ment,  that  is,  if  53i=i  ni ,i-i  ~  N)  or,  since  the  Wi>t_i(<)  are  not  integers,  when 

(4.10)  N  -  ±  nj.i-i  =  ±  =  Ne~”»  =  N  -  ~  ”  ,  <  y 

In  the  stochastic  case,  if  we  write 

s  —  1 

(4.11)  72  =  72(r?00,  n0i,  •  •  •  ,  n0s_ ij  ni.;  £),  »i-  =  N  —  £  n0i, 

i  =  0 

as  the  probability  that  there  are  at  time  t,  n0o,  n0i,  •  •  •  ,  n0s-i  phages  with  0,  •  •  •  , 
s  —  1  attachments,  respectively,  at  other  than  tail  positions,  then,  assuming 
for  x (t)  the  deterministic  solution  (2.11),  we  have  that 

(4.12)  f  = 

8—1 

^ '  ot(s  i  l)(woi  “l-  1 ) xR (n00,  ,  n0»  — {—  1,  Woi+i  1,  *  *  *  ,  nos— i ,  ni. ,  2) 

t=0 

8—1 

+  53  «(«oi  +  l)xR(n00,  •  •  •  ,  n0;  +  1,  n0i+i,  •  •  •  ,  n0s_i;  n\.  —  1;  t) 

i  =  0 
8—  1 

—  53  «(«  —  i)n<sixR{nm,  •  •  •  ,  n0»_i;  ni.;  /)• 

i  =  0 


540 


FIFTH  BERKELEY  SYMPOSIUM:  GANI 


The  generating  function  for  these  probabilities  w0o ,  •  •  •  ,  ;  v ;  f)  satisfies 

the  partial  differential  equation 

(4.13)  “7  =  E  ax{(s  —  i  —  1)m0i+i  +  v  —  (s  —  i)u o*} 

at  i=o 


We  may  write  in  the  usual  way  that 

dt  d\[/  dv  duoo 


(4.14) 


du0s — x 


1  0  0  a:r[(s  —  l)w0i  +  v  —  suoo] 

so  that  if  U'  =  ( Uw ,  •  •  •  ,  Wo«-i),  then 

s  -(s-l) 


ax[v  —  swo(s-i)] 


(4.15)  ^  =  otx 


-1)  — («-  2) 


1 


U  —  ax  lv 


=  ax(UJ  -  yl). 

Treating  y  as  a  constant,  the  solution  of  this  is  seen  to  be  of  the  form 

(4.16)  e-«p(OL(lJ  —  yL-1l)  =  constant, 

and  it  follows  that 


(4.17)  ^(woo,  •  •  •  ,  t)  =  *[e  “p(0L( U  —  yL  *l),y] 

subject  to  the  condition  that  ^(woo,  •  •  •  ,  u0s]  y;  0)  = 

Thus, 

^(U  —  yL-1l,y)  =  Uoo, 


so  that,  since  L_1l  =  1  for  the  matrix  L  in  (4.15),  we  have 

(4.18)  *[e-«p(t)L(tJ  _  yL-1!),  y]  =  {[e-“p«L(U  -  yL-1l)]0  +  y}*, 


where  [e-ap(‘)L(U  —  yL_1l)]0  indicates  the  Oth  element  of  the  column  vector. 

It  may  readily  be  shown  after  some  matrix  calculations  that  the  probability 
generating  function  is  of  the  form 


(4.19) 


iA(woo,  >  ^os— ij  y  f  t) 


£ 

i  =0 


where  the  &,•(<)  are  of  the  form 

b0(t )  =  e~,ap(-t)  =  {e~pt  +  s(s  —  —  e_#x0}~*> 

l  bi(t)  =  (S  7  1)  e_(8_l)ap(0(l  —  e-°p(*))‘,  i  =  1,  •  •  •  ,  s  —  1. 

These  probabilities  add  up  to 

(4.21)  *E  T  *)  e-(8_i)ap(l  -  e~apy  =  e~ap(t). 

Thus,  for  phage  virus,  the  probability  of  loss  of  infectivity  is  given  by 

r  in  ( i  —  c  ^ 

(4.22)  P{n m  =  0,  •  •  •  ,  n0._,  =  0}  =  (1  -  e— 0*  -  (  ,  -  me~''  J  ' 
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5.  Conclusion 

Further  investigations  of  more  realistic  models  and  their  verification  in  the 
laboratory  would  be  of  interest.  In  the  case  of  the  influenza  virus,  for  example, 
it  is  known  that  antibodies  may  bend  over  to  attach  both  their  ends  to  emplace¬ 
ments  on  the  same  virus  particle.  It  is  also  possible  for  one  end  of  an  antibody 
to  be  attached  to  an  emplacement  on  one  virus,  while  its  other  end  is  attached  to 
a  second  virus  particle;  conglomerations  of  viruses  and  antibodies  can  thus  be 
formed.  Clearly,  the  geometry  of  such  models  becomes  more  complicated  than 
that  we  have  outlined  earlier. 

It  may  still  be  possible,  however,  to  construct  simplified  models  for  them, 
and  to  draw  probabilistic  conclusions  from  these.  Similarly,  for  bacteriophage, 
a  model  could  be  constructed  in  which  antibodies  have  one  end  attached  to  the 
tail  of  one  phage  particle  while  the  other  is  attached  to  a  second  phage  tail* 
Such  a  model  is  not  too  intractable,  and  it  is  hoped  to  present  results  relating 
to  it  in  some  work  at  present  in  preparation. 


REFERENCES 

[1]  M.  S.  Bartlett,  “Some  evolutionary  stochastic  processes,”  J.  Roy.  Statist.  Soc.  Ser.  B, 
Vol.  II  (1949),  pp.  211-229. 

[2]  J.  Gani,  “Stochastic  phage  attachment  to  bacteria,”  Biometrics,  Vol.  21  (1965),  pp. 
134-139. 

[3]  - ,  “Stochastic  models  for  bacteriophage,”  J.  Appl.  Probability,  Vol.  2  (1965),  pp. 

225-269. 

[4]  E.  N.  Gilbert,  “The  probability  of  covering  a  sphere  with  N  circular  caps,”  Biometrika , 
Vol.  52  (1965),  pp.  323-330. 

[5]  P.  A.  P.  Moran  and  F.  de  St.  Groth,  “Random  circles  on  a  sphere,”  Biometrika,  Vol.  49 
(1962),  pp.  384-396. 

[6]  W.  L.  Stevens,  “Solution  to  a  geometrical  problem  in  probability,”  Ann.  Eugen.,  Vol.  9 
(1939),  pp.  315-320. 

[7]  D.  Yassky,  “A  model  for  the  kinetics  of  phage  attachment  to  bacteria  in  suspension,” 
Biometrics,  Vol.  18  (1962),  pp.  185-191. 


